Infrared structure of e + e —> 3 jets at NNLO: QED-type contributions 
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The NNLO QCD corrections to the e + e~ — » 3 jets can be decomposed according to their colour factors. Out of 
the seven colour factors, three are of QED-type: 1/N 2 , N F /N and N 2 F . We use the antenna subtraction method 
to compute these contributions, providing complete expressions for the subtraction terms in Nf/N and N F . 



1. Introduction 

Electron-positron annihilation into three jets 
provides one of the most precise experimental 
tests of QCD dynamics. The extraction of the 
strong coupling constant a s from data on this ob- 
servable 1 is currently limited 2 by the uncer- 
tainty on the available theoretical next-to-leading 
order (NLO) calculations. To improve upon this 
situation, it is mandatory to compute next-to- 
next-to-leading order (NNLO) QCD corrections 
to e + e~ — > 3 jets. 

This calculation contains three ingredients [3]: 
matrix elements for 7* — > 5 partons at tree-level, 
for 7* — > 4 partons at one-loop and 7* — > 3 par- 
tons at two loops. The four-parton and five- 
parton contributions contain infrared singulari- 
ties due to real radiation of massless final state 
partons. To extract these singularities, which 
cancel with explicit infrared singularities in the 
virtual corrections for any infrared-safe observ- 
able, one needs to define a subtraction scheme at 
NNLO. We proposed a formulation of the antenna 
subtraction method gE] to NNLO in 

2. Antenna subtraction 

In the antenna subtraction method, antenna 
functions describe the colour-ordered radiation of 
unresolved partons between a pair of hard (radi- 
ator) partons. All antenna functions at NLO and 
NNLO can be derived systematically from phys- 
ical matrix elements, as shown in |6I7) . They 



can be integrated over the factorized antenna 
phase space jS] using loop integral reduction tech- 
niques extended to phase space integrals |9I10| . 
and then combined with virtual corrections to 
partonic processes with lower multiplicity. 

3. Colour structure of e + e — > 3 jets 

Decomposing the three-jet cross section (and 
related event shape observables) at NNLO ac- 
cording to the QCD colour structures one finds 
seven colour factors: A 2 , A , 1/A 2 , N F N, 
Np/N, Np and Np }1 . The last colour factor 
Ap j7 arises from the interference of amplitudes 
with two independent quark lines coupling to the 
exchanged vector boson. This colour factor is in- 
frared finite and numerically small in the three- 
parton channel and contributes only a negli- 
gible amount to four-jet observables at NLO |12| . 
It can therefore be neglected safely. 

Among the remaining six colour factors, three 
do not receive contributions from diagrams with 
gluon self-coupling, and are therefore of QED- 
type: 1/A 2 , N F /N and Np. Tn we de- 
scribed the calculation of NNLO corrections in 
the 1/A 2 colour factor using the antenna sub- 
traction method. In the present paper, we de- 
scribe the calculation of NNLO corrections in the 
Np/N and Np colour factors, thus completing 
the corrections to all QED-type colour factors. 
Our notation for partonic matrix elements, an- 
tenna functions and subtraction terms follows the 
notation introduced in [Sj. 
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4. Subtraction in the Np/N colour factor 

The Np/N colour factor receives contributions from five-parton tree-level 7* — > qqq'q'g, four-parton 
one- loop 7* — > qqq'q' and 7* — > gggg at subleading colour as well as tree- level two- loop 7* — > qqg. The 
gluon emissions are all photon-like. 

4.1. Five-parton contribution 

The NNLO radiation term appropriate for the three jet final state is given by 

d(J NNLO.N F /N = N 5 ^d$ 5 (pi,...,p 5 ;gr) B°' c (l q ,5 g ,2 q ;3 q ,,4g,) 

+B 5 M (1 9 , 2 9 -; 3 q ,, 5 g , 4 r ) - 2 £? 5 °' e (l 9 , 2 ff ; 3,,,4 r ; 5„)1 jfV, . . . ,p B ) 



= N 5 ^d$ 5 

(i,j')e(3,4) 



(l„2g; v,5 fl , j>) - 2B°- e (l q ,2q;i q ,, jq<;5 g ) J^'ipi, ■ ■ ■ ,p 5 ) , (1) 

where the symmctrisation over the momenta of the secondary quark-antiquark pair exploits the fact that 
the jet algorithm does not distinguish quarks and antiquarks. This symmctrisation reduces the number 
of non-vanishing unresolved limits considerably, since the interference term in B®' e is odd under this 
interchange. The subtraction term reads: 

o Np 

d(T NNLO.N F /N = N ^^~ d$ 5(Pl, • • • ,P5 



6(3,4) I 



-A§(l„5 fl ,2 ff )B2((15),,(25) ff ,v,i r )J r i 4) (^,p 2 B,Pi,Pj) 

-A% (i q , , 5 g , jg, )Bl(l q ,2q,(i5) ql , (j5) ql ) f 3 4) (j>l,P2,Pi5,Pj5) 

- \ {^3° (1 9 , V , 39 ) Al{ilj)g, W) g , 5 fl) 2,-) Ji 4) (py , pa, p£, P5) + (1 <- 2) } 

-^ 4 (l g , V,jV,2 ff ) - \ {E°(lgJg,Jg,)Al((lj) q M g ,%) + ^ " ^ ^ 

x^3((W) 9 , 5 fl , (2ij) q ) 4 3) (Pw,P^j,P5) 
+ \^El{l q ,ig,dg,)Al{(lJ)g,b g ,2q) A°(((tf)5),,(tf) fl) (25)-) jf (^ 5 ,pT 5 ,^) + (1 - 2) I 

xA°((l^) 9 ,(5^) g ,2 (? )jf ) (^,p 2 ,^) + (1^2)|| (2) 

4.2. Four-parton contribution 

The four parton contribution to the Nf/N colour factor reads: 

d<J NNLO.N F /N = N 4^j- (f^) d* 4 (pi,...,P4;«) |- £ E (^''(^.V.JV. 2 ?) 

L (»,j)€(3,4) 
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+2Cl J (l q ,z q , Jq ,2 q ) + A^(l q , ig , Jg ,2 q )) } Jf>(p 



The expression is symmetrised over the momenta (3) and (4) to remove terms which are antisymmet 
under charge conjugation, and can not be accounted for properly by the quark-gluon antenna funct ; "' 
The corresponding subtraction term is: 



(3) 

;ric 
ions. 



+ \ [£l(s 13 ) + £l{s u ) + El{s 2Z ) + Zl{ S2i )] A 4 (l q ,3 g A g ,2 q )4 4) (Pu.-.,P4) 

-^UE° s (l q ,3 q ,A?)[Al((U) q ,m 
+A 3 (s ( ~ )2 ) E°(l q , 3 q ,,4,,) A° 3 ((13) q , (34) s , 2 q ) 

El(l q ,3 q ,,4 q ,) +A Q 3 (s 34 )E° 3 (l q ,3 q ,, 4 q ,)] A° 3 ({13) q , (34) g ,2 q )\ jf (pTa.pSi.Pa) + (1 - 2) 



~ E (^3(l g »V2 9 )i5((l»),,i fl ,(2i) ff ) + 

(i,j)e(3,4) 

1 



+ - (£°( Sli ) + £$( Slj ) + £° 3 (s 2t ) + £° 3 (s 2j )) A° 3 (l q , i g , 2 q ) A° 3 ((li) q ,j g , (2i)-) 



+2b , F log^-A 3 (l q ,i g ,2 q )A° 3 ((li) q ,j g ,(2i) ))j< 3) (pTufcuPj) 
sin H H 



(4) 



The terms in the second and third line of this equation, as well as in the second-last line are obtained by 
integrating terms from the five-parton subtraction term @. All remaining terms are integrated to yield 
three-parton contributions. 

4.3. Three-parton contribution 

The three parton contribution to the Np /N colour factor consists of the three-parton virtual two-loop 
correction and the integrated five-parton tree-level and four-parton one-loop subtraction terms, which 
read 



da 



NNLO,N F /N + & a NNLO,N F /N ~ 

1 =n , . 1 = 

2 



vs,i 



N,- 



x< - 



B° 4 (s 12 ) - i^(a 13 ) -\s° 4 ( S23 ) - ^AtM (£l( Sl3 ) + £° 3 ( S23 )) 
>i 2 ) + l^is) + \el{s 2 s) A°(l„3 fl ,2,-) - l - (£%s 13 ) +£° 3 (s 23 )) A*(l„3 flJ 2 a -) 
A%s 12 )Al(l q ,3 g ,2 q ) - b -^A%s 12 ) (( Sl2 )- e - ( Sl23 )- £ ) A°(l 9 ,3 s ,2 9 -)|da 3 



(5) 



Taking the infrared pole part of this expression, we obtain cancellation of all infrared poles in this 
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channel: 



Voles (da s NNLO NF/N ^j + Voles (dCT^ O AWAr ) Voh * 
5. Subtraction in the Np colour factor 



(da 



V,2 

NNLO,N F /N 



(6) 



The Np colour factor receives contributions only from the four-parton one-loop process 7* 
and from the three-parton two- loop process 7* — > qqg. 

5.1. Four-parton contribution 

The four-parton one-loop contribution to this colour factor is 



qqq'q' 



da 



V,l 

NNLO,Nl 



N±N F (g) d*4(pi,...,P4;g)B 4 1,c (l„3,,,4 r ,2 ? ) ji 4) ( Pl ,..., P4 ). 



(7) 



This contribution is free of explicit infrared poles (as can be inferred from the absence of a five-parton 
contribution to this colour structure). 

The subtraction term appropriate to this contribution is 



A(jV n S nlo,ni = n * n f (^) d$ 4 (Pi,...,P4;g)5 



£ 3 1 (l„3,-,4 r ) + 2VF log ^£ 3 °(1 9 , 3^,4,0 A°((13) gI (34) g) 2 ? ) 

S134 



+^(l 9 ,3,,4 r )i5((13),,(34) fl ,2 9 ) 



^i 3) (Pl3,P34,P2) 



E\{2 9 ,\^) + 2b 0>F log^£; 3 (2 g -,3 gS 4 9 -0 A§(1„ (34),, (23) ? ) 

S234 



+E° 3 (2 q , 3 q >, 4,-) i*(l„ (34),, (23).-) 



4 3) (^l'-P34,P23) 



(8) 



Although S3 and A3 contain explicit infrared poles, these cancel in their sum, as can be seen from (5.16) 



and (6.32) of 0. da 



VS,1 

NNLO,Nl 



is therefore free of explicit infrared poles. 



5.2. Three-parton contribution 

The three parton contribution to the Np colour factor consists of the three-parton virtual two-loop 
correction and the integrated four-parton one-loop subtraction term, which reads 



da 



vs,\ 



Ni 



NNLO,N% ~ r 2 

4 x ( Sl3 ) + 4 X (*23)) A° 3 (l q , 3 g , 2 q ) + (£ 3 °( Sl3 ) + £ ( i(s 23 )) Al(l q , 3 g , 2 q ) 
+ b -^f [e$(s 13 ) (( Sl3 )- £ - ( Sl23 r) +£$(s2s) (( S23 )- £ - ( S i 23 )- £ )] A°(l 9 ,3„2 ? ) 



d(7 3 



(9) 



Taking the infrared pole part of this expression, we obtain cancellation of all infrared poles in this 
channel: 



Voles 



(a vs,i 

{ d<T NNL 



NNLO,Ni 



Voles [ da 



( da N 2 NLO,Nl) = °- 



(10) 
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6. Numerical implementation 

Starting from the program EERAD2 4 , which 
computes the four-jet production at NLO, we 
implemented the NNLO antenna subtraction 
method for all QED-type colour factor contribu- 
tions to e + e~ — » 3j. EERAD2 already contains the 
five-parton and four-parton matrix elements rel- 
evant here, as well as the NLO-type subtraction 
terms. 

The implementation contains three channels, 
classified by their partonic multiplicity: (a) in 
the five-parton channel, we integrate A^nnlo ~ 
d&NNLO' ( D ) m t ne four-parton channel, we inte- 
grate dcr^ io — dcrj^^ ; (c) in the three-parton 
channel, we integrate dcrJy^ LO + dcr^ NLO + 

V S 1 

^nnlo- ^he numer i ca l integration over these 
channels is carried out by Monte Carlo methods. 

By construction, the integrands in the four- 
parton and three-parton channel are free of ex- 
plicit infrared poles. In the five-parton and four- 
parton channel, we tested the proper implementa- 
tion of the subtraction by various checks. To ver- 
ify the local convergence of the subtraction terms 
and to verify that cancellations among individ- 
ual contributions to the subtraction terms take 
place as expected, we examined the behaviour 
of matrix element and subtraction terms along 
various trajectories approaching the unresolved 
limits. Moreover, we checked the correctness of 
the subtraction by introducing a lower cut (slic- 
ing parameter) on the phase space variables, and 
observing that our results are independent of this 
cut (provided it is chosen small enough). This 
behaviour indicates that the subtraction terms 
ensure that the contribution of potentially sin- 
gular regions of the final state phase space does 
not contribute to the numerical integrals, but is 
accounted for analytically. Finally, distributions 
in double and triple invariants of the five-parton 
or four-parton phase space illustrate the proper 
onset of the subtraction terms towards the single 
and double unresolved edges of phase space. 

7. Outlook 

In this paper, we described the calculation of 
NNLO corrections to e + e~ — > 3 jets in all QED- 
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type colour factors: l/N 2 , N F /N and Np, us- 
ing the antenna subtraction method [H] ■ We doc- 
umented explicit expressions for the subtraction 
terms for N F /N and Np (l/N 2 was documented 
in [d] already) and described their implementa- 
tion in a numerical program computing three-jet 
cross sections and related event shape observ- 
ables. Work on the remaining colour factors is 
ongoing. 
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